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EVIDENCE FOR ACCRETION RATE CHANGE DURING TYPE I X-RAY BURSTS 
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ABSTRACT 

The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts involves subtraction 
of the pre-burst emission as background. This approach implicitly assumes that the persistent flux remains 
constant throughout the burst. We reanalyzed 332 photospheric radius expansion bursts observed from 40 
sources by the Rossi X-ray Timing Explorer, introducing a multiplicative factor f a to the persistent emission 
contribution in our spectral fits. We found that for the majority of spectra the best-fit value of f a is significantly 
greater than 1, suggesting that the persistent emission typically increases during a burst. Elevated f a values 
were not found solely during the radius expansion interval of the burst, but were also measured in the cooling 
tail. The modified model results in a lower average value of the x 2 fit statistic, indicating superior spectral 
fits, but not yet to the level of formal statistical consistency for all the spectra. We interpret the elevated f a 
values as an increase of the mass accretion rate onto the neutron star during the burst, likely arising from the 
effects of Poynting-Robertson drag on the disk material. We measured an inverse correlation of f a with the 
persistent flux, consistent with theoretical models of the disc response. We suggest that this modified approach 
may provide more accurate burst spectral parameters, as well as offering a probe of the accretion disk structure. 



1. INTRODUCTION 

Thermonuclear (type I) X-ray bursts arise from the unsta- 
ble ignition of accreted H/He near the surface of an accreting 
neutron star in a low-mass X-ray binary (LMXB) (e.g., Fuji- 
moto et al. 1981; Strohmayer & Bildsten 2006). Gas accreted 
from a low-mass stellar companion accumulates on the sur- 
face of the neutron star, where it is compressed and heated 
hydrostatically. When the temperature and pressure are high 
enough a thermonuclear explosion is triggered. These events 
are observed as a sudden increase in X-ray luminosity to many 
times the persistent level (e.g., Lewin et al. 1993; Strohmayer 
& Bildsten 2006). Type I bursts have been detected from over 
100 sources in our Galaxy 4 . Typical bursts exhibit rise times 
of 1 to 10 seconds, durations of a few tens of seconds to a few 
minutes and have total energy outputs of 10 39 - 10 40 erg. 

Analyses of type I bursts typically make a number of im- 
plicit assumptions, namely: 

• The total source spectrum consists of two additive com- 
ponents: one (the "burst component") arising from nu- 
clear burning and the other ("persistent emission") aris- 
ing from accretion. 

• The burst component has the same spectral shape for all 
bursts from all sources. 

• That shape is a blackbody, with variable temperature 
and normalization. 

• The persistent emission's spectral shape does not 
change during a burst and is identical to its pre-burst 
shape. 
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• The persistent emission's intensity does not change dur- 
ing a burst and is identical to its pre-burst intensity. 

Making the last two of these assumptions allows the sub- 
traction of the pre-burst emission of the neutron star from 
the burst spectra as background (e.g., van Paradijs & Lewin 
1986; Lewin et al. 1993; Kuulkers et al. 2003; Galloway et al. 
2008, hereafter G08). This approach ( hereafter referred to 
as the "standard approach") implicitly assumes that the accre- 
tion rate remains constant throughout a burst, but it is not ob- 
vious that this assumption is reasonable. For instance, when 
the flux reaches the Eddington limit, in the so-called photo- 
spheric radius expansion (PRE) bursts, one would naively ex- 
pect, assuming isotropic emission, accretion to cease entirely 
because the outward radiation force exceeds the gravitational 
force. Lamb & Miller (1995) argue that accretion ought to be 
shut off if the luminosity exceeds Eddington anywhere in the 
accretion flow, not necessarily just at the stellar surface. The 
effect of radiation halting accretion flow has apparently been 
observed for a number of the very brightest and most vigor- 
ous bursts (e.g., in't Zand & Weinberg 2010), during which, 
for a few seconds, no X-ray flux (beyond instrumental back- 
ground) is observed from the star — though this may also be 
due to the atmosphere expanding so far that it obscures the 
emitting regions of the accretion disc. On the other hand, in- 
creased luminosity during a burst might enhance the accretion 
rate, via Poynting-Robertson drag acting upon the accretion 
disc (Walker & Meszaros 1989; Miller & Lamb 1996; Walker 
1992, hereafter W92). At luminosities greater than O.OlLEdd, 
radiation forces have more of an effect on the accretion flow 
than general relativistic effects (Miller & Lamb 1993). Since 
there are questions about the isotropy of the burst emission 
(e.g., Boutloukos et al. 2010), it is not obvious which of out- 
wards pressure and radiation drag will dominate. It is there- 
fore of importance to determine if a varying accretion rate i s 
detectable during a PRE burst, and to quantify any variation 
that is detected. In this paper we attempt to measure a change 
in the accretion rate by performing spectral fits where the pre- 
burst (persistent) emission is allowed to vary during a burst, 
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but holding all the other assumptions of the standard approach 
fixed. 

Some of the implicit assumptions of the standard analysis 
approach have been tested in previous studies: van Paradijs & 
Lewin (1986) pointed out that, if the total spectrum and per- 
sistent spectrum are both treated as blackbodies, subtracting 
the former from the latter will not leave a net burst spectrum 
that can be fit with a blackbody. This idea was followed up 
by Kuulkers et al. (2002) in a study of the rapidly accreting 
GX 17+2. They did not find that accounting for this effect 
improved the spectral fits and concluded that the persistent 
emission does not originate from the same location as the 
burst emission on that neutron star. Muno et al. (2000) and 
Strohmayer & Brown (2002) allowed for the accretion to shut 
off entirely, by subtracting the instrumental background only, 
but did not find that the spectral fits were improved by do- 
ing so. Recently in't Zand et al. (2013) studied a PRE burst 
from SAX J1808. 4-3658 using combined Chandra and RXTE 
data and found that an observed excess of photons at both low 
and high energies can be well described by allowing a 20-fold 
increase of the pre-burst persistent emission. It also may be 
that the persistent emission is composed of separate contribu- 
tions arising from different sites. These include a boundary 
layer at the inner edge of the accretion disk (e.g., Kuulkers 
et al. 2002), the inner regions of the accretion disk proper 
(e.g., Christian & Swank 1997; Cackett et al. 2010), emission 
from the neutron star itself or its photosphere (van Paradij s 
& Lewin 1986), and Compton scattering in an accretion disk 
corona (White & Holt 1982). However, disentangling these 
contributions is likely to be difficult, because they are proba- 
bly correlated and spectrally indistinct. 

Deviations of the burst component of the spectrum from 
a black-body spectrum could also be present. Such spectral 
changes have been theoretically predicted at both the high en- 
ergy (e.g., London et al. 1984, 1986), and low energy (e.g., 
Madej 1991; Madej et al. 2004) ends of the X-ray spectrum. 
However, the literature is divided as to whether these are ac- 
tually present in observations. Excess photons at high energy 
have been reported in bursts from 4U 2129+1 1 (van Paradijs 
et al. 1990) and GX 17+2 (Kuulkers et al. 2002). On the 
other hand, pure black-bodies have been found to give gen- 
erally good results up to the present time (e.g., Giiver et al. 
2012) and some authors argue that they are consistent with 
black-bodies to extremely high confidence (Boutloukos et al. 
2010). One might further divide the burst emission into con- 
tributions from a continuum, and discrete spectral features 
superimposed upon it such as emission lines and absorption 
edges. Continuum changes are likely to be present throughout 
all stages of a burst (e.g., Suleimanov et al. 2011). Changes 
in the spectral features are thought to be largely confined 
to the Eddington-limited radius expansion period as these 
are thought to be due to ashes from nuclear burning being 
mixed into the expanding envelope (Weinberg et al. 2006), 
and have been detected in the so-called superexpansion bursts 
(in't Zand & Weinberg 2010). However, the non-detection 
of spectral features in a PRE burst from SAX J1808.4-3658 
by in't Zand et al. (2013) suggests that such features may be 
too weak to be detected by currently available instruments: 
this source is the brightest PRE burster (see Table 4) and the 
burst in question was observed by two X-ray observatories. 
Galloway et al. (2010b) previously used Chandra spectra to 
search for spectral features in PRE bursts from 4U 1728-34, 
without success. Either the sought-after features are too weak 
to be detected, or they are not present in every burst. 



If the accretion rate varies for any bursts, we expect it to be 
for the PRE bursts, although there is obviously a possibility 
that outward radiation pressure and radiation drag will par- 
tially negate each other; this makes them the most stringent 
tests of the idea that accretion rate might increase. They are 
the most luminous bursts for any given source and therefore 
should produce the largest radiation forces. They will also 
have the best signal-to-noise of bursts from any given source. 
These bursts are of interest because they can be used to probe 
the structure of LMXB systems. Since they reach the Edding- 
ton luminosity of the neutron star, the atmosphere is no longer 
bound to the surface of the star and expands. The luminosity 
is thought to remain within a few percent of the local Edding- 
ton luminosity throughout the radius expansion (e.g., Hanawa 
& Sugimoto 1982; Ebisuzaki et al. 1983; Titarchuk 1994b). 
Thus, radius expansion bursts can in principle be used to mea- 
sure the surface gravitational redshift of the neutron star and 
thereby its compactness (e.g., Damen et al. 1990; Ozel et al. 
2009; Giiver et al. 2010; Steiner et al. 2010), giving insights 
into the neutron star equation of state (Lattimer & Prakash 
2007). Radius expansion bursts can also be used to determine 
the distance to the bursting star (Basinska et al. 1984), mak- 
ing them potentially useful as standard candles (e.g., Kuulkers 
et al. 2003). 

The structure of this paper is as follows. In §2 we describe 
the data and its collection. In §3 we develop and implement 
a modified spectral analysis method, and apply it to a single 
burst to demonstrate its feasibility. In §4 we apply the method 
to all the PRE bursts in our sample to build a statistical picture 
of the phenomenon. In §5 we investigate the effect of allowing 
non-thermal burst emission of fixed shape. In §6 we compare 
our results to theoretical predictions made by Walker (1992). 
In §7 we discuss our results and place them in the context of 
previous work. 

2. DATA COLLECTION AND REDUCTION 

We used observational data from the Rossi X-Ray Tim- 
ing Explorer (RXTE), publicly available through the High- 
Energy Astrophysics Science Archive Research Centre 
(HEAS ARC), dating from shortly after the satellite's launch 
on December 30, 1995 to the end of the RXTE mission on Jan- 
uary 3, 2012. RXTE carries three instruments for detecting 
X-rays. The All-Sky Monitoring camera (ASM) consists of 
three Scanning Shadow Cameras sensitive to photons with en- 
ergies between 1.5 and 12 keV with a field of view of approx- 
imately two degrees, and performed 90s step-stare observa- 
tions of most of the sky every 96 minutes (Levine et al. 1996). 
The Proportional Counter Array (PCA) consists of five PCUs 
sensitive to photons with energies between 2 and 60 keV and 
has a field of view of approximately one degree (Jahoda et al. 
2006). RXTE also carries the High Energy X-ray Timing Ex- 
periment (HEXTE), a collection of scintillation counters with 
a one degree field of view (Rothschild et al. 1998), but we do 
not use data from HEXTE in this work. 

We took as the basis for our sample the G08 catalog of 
bursts. The G08 sample identifies PRE bursts according to 
the criteria set out in G08 (their section 2.3). Briefly, these 
criteria define a PRE burst as one that i) reaches a local max- 
imum blackbody normalization at or near the moment of 
peak flux, ii) has declining values of after this time, and 
iii) attains its lowest blackbody temperature at the maximum 
A"bb. At the time of publication of G08, the catalog contained 

5 See http://heasarc.gsfc.nasa.gov 
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254 PRE bursts; since then, a further 118 PRE bursts have 
been observed, giving a total of 372 6 . 

Unless otherwise stated, the data analysis procedures are as 
in G08. For those bursts for which suitable datamodes were 
available, time-resolved spectra in the range 2-60 keV cover- 
ing the burst duration were extracted on intervals beginning at 
0.25 s during the burst rise and peak. The bin size was gradu- 
ally increased into the burst tail to maintain roughly the same 
signal-to-noise level. The RXTE Proportional Counter Units 
(PCUs) are subject to a short (« 10 ps) interval of inactivity 
following the detection of each X-ray photon. This "dead- 
time" reduces the detected count rate below what is incident 
on the detector (by approximately 3% for an incident rate of 
400 count s" 1 PCU" 1 ). We calculated for each measured spec- 
trum an effective exposure, taking into account the fraction of 
events lost during deadtime, following the approach recom- 
mended by the instrument team 7 . Contributions to the dead 
time fraction arise from coincidence and particle events as 
well as source and background photons. 

We re-fit the spectra over the energy range 2.5-20 keV us- 
ing the revised PCA response matrices, vll.7 8 and adopted 
the recommended systematic error of 0.5%. The fitting was 
undertaken using XSpec version 12. In order to accommo- 
date spectral bins with low count rates, we adopted Churazov 
weighting. 

We modelled the effects of interstellar absorption, using a 
multiplicative model component (wabs in XSpec), with the 
column density A^h frozen. The A^ values are drawn from the 
literature, preferentially from studies of neutron stars using in- 
struments sensitive at lower X-ray energies than RXTE. They 
are listed in Appendix A. In the original analysis carried out 
by G08, the neutral absorption was determined separately for 
each burst, from the mean value obtained for spectral fits car- 
ried out with the A^h value free to vary. This has a negligible 
effect on the burst flux, but can introduce spurious burst-to- 
burst variations in the blackbody normalization. 

3. METHOD 

Our revised analysis is to fit the burst spectra with a two- 
component model consisting of a black-body, representing the 
burst emission, and a model of the pre-burst persistent emis- 
sion, representing the emission due to accretion of material 
onto the neutron star, with a prefactor f a left free in the fits. As 
we cannot distinguish contributions to the persistent emission 
arising from different locations in the neutron star system, we 
simply assume the persistent emission is indivisible and re- 
sults entirely from accretion onto the neutron star. As PRE 
bursts are likely to be the most stringent tests of changes in 
accretion onto the star, we select these events for our analysis. 

It is also possible that the burst emission is non-Planckian, 
but fitting a non-black-body spectrum requires the existence 
of theoretical models that can describe the data. The most re- 
cent model atmospheres of bursting neutron stars are those of 
Suleimanov et al. (201 1), but even these have been only par- 
tially successful (Zamfir et al. 2012), and are not intended for 
radius expansion spectra. In the absence of models that are 
demonstrably better than black-bodies, and no consensus that 
black-body fits really are unsuitable, we keep the assumption 
of thermal burst emission for the majority of our analysis. We 

6 see burst.sci.monash.edu/minbai' 

7 see http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pca_deadtime.html 

8 see http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/pcarmf- 
11.7 



discuss this issue further in §5. Similarly, we do not draw a 
distinction between deviations from a blackbody due to con- 
tinuum or features in this paper, since spectral features are ei- 
ther not present in every burst or are too weak to be detected. 

In order to demonstrate our approach, we select a burst 
from the well-studied PRE burster 4U 1636-536 (Swank 
et al. 1976). This source is an attractive candidate for sev- 
eral reasons. It is a prolific burster, with 75 PRE bursts 
recorded by RXTE. At a distance of approximately 6 kpc 
(G08) it is relatively nearby. Its average peak burst flux of 
(69 ± 6) x 10~ 9 erg cm" 2 s" 1 is among the brightest sources 
and ensures good signal-to-noise. The hydrogen column den- 
sity of 0.25 x 10 22 cm" 2 (Asai et al. 2000) is low compared to 
most other sources in the catalog, minimizing the absorption 
corrections that have to be performed on the spectra. There 
are no other known LMXBs in the same field of view as 
4U 1636-536, so confusion with the persistent emission of 
other sources is not an issue (see §4.1 for further discussion 
of this problem). Finally, 4U 1636-536 accretes mixed H/He 
(Galloway et al. 2006) and is therefore representative of the 
majority of neutron stars in our sample (G08). 

The burst we chose was detected on Jun 15, 2001 (burst ID 
34 in G08). The data for this burst consists of 122 spectra 
recorded by RXTE with an integration time of typically 0.25s. 
Of these spectra, 37 are after the beginning of the burst, which 
we take to be the first time that the flux exceeds 25% of the 
maximum burst flux. Recording continued up to 176.75s after 
the burst start, but the flux had declined to the pre-burst level 
by 24.75s. Two of RXTE's PCUs, numbers 3 and 4, were 
active for this burst. 

3.1. Characterizing the persistent emission 

We adopted the integrated X-ray flux for a 16-second in- 
terval prior to the start of each burst as the persistent emis- 
sion. This spectrum also includes a time-dependent contribu- 
tion from instrumental (non-source) background; to estimate 
this contribution we used the full-mission, "bright" source 
(> 40 counts" 1 ) models released 2006 August 6 with the 
pcabackest tool. Subsequent model fits to each persistent 
(and burst) spectrum used the corresponding model spectrum 
estimated for that burst as background. We then generated a 
model for the persistent emission consisting of a blackbody 
plus a power law, both corrected for interstellar absorption. 
The fits were performed with XSPEC (Arnaud 1996; Dorman 
& Arnaud 2001) using wabs* (bbodyrad+powerlaw) as 
the model. The hydrogen column density was kept fixed at 
0.25 x 10 22 cm" 2 , the value determined in Asai et al. (2000). 

For the persistent emission model we obtained a black- 
body temperature of (1.9_q[) keV and a normalization of 
(6.0_[q) km 2 /(10 kpc) 2 . The powerlaw component had index 
(+3.0!^) and normalization (2.2!°j)keV- 1 cm- 2 s" 1 at 1 keV. 
Errors given are the 1 a confidence level. The reduced % 2 for 
this fit was 0.699 for 21 degrees of freedom, indicating that 
the model adequately describes the data. Figure 1 shows the 
fit to the data and the residuals. 

3.2. Modelling the burst spectrum 

We initially fit the burst spectra with the standard ap- 
proach: a blackbody spectrum, corrected for interstellar ab- 
sorption using the same hydrogen column density as above, 
ie. wabs*bbodyrad in XSPEC. We subtracted the mea- 
sured pre-burst emission, which includes a component that 
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FIG. 1. — Measured persistent spectrum for the Jun 15, 2001 burst from 
4U 1636—536 and the fitted spectral model (top panel), and residuals to the 
data (bottom panel). The reduced \ 2 is 0.699, indicating a good fit to the 
data. 



does not arise from the source (the instrumental background), 
and fit the resulting burst spectrum. This is the same imple- 
mentation of the standard approach as in G08. These initial 
fits provide standard approach fits to which we can compare 
ours, as well as sensible spectral parameters to seed the vari- 
able persistent emission fits. 

Now we introduce the dimensionless quantity f a , the per- 
sistent emission multiplicative factor, as a third variable pa- 
rameter. The burst spectra are fitted again, replacing the 
recorded pre-burst emission subtraction with just the instru- 
mental background, with the model 

S(E)=A(E) x B(E;T bb ,K bb )+f a x P(E) - b(E) msU (1) 

where S(E) is the fitted spectrum as a function of energy E, 
A is the absorption correction, B is the blackbody function 
with temperature T bb and normalization K bb , P is the persistent 
model described in §3.1, and Z?j nst is the instrumental back- 
ground. Note that the persistent model already contains an 
absorption factor. The parameter f a is allowed to vary be- 
tween -100 and 100. This allows us to track changes in the 
accretion rate through the factor f a . Because of the 0.25 s 
integration time of the burst spectra, our measured peak f a 
values may be slight underestimates. 

Figure 2 shows the time evolution of burst bolometric flux 
(top panel), the f a factor (middle panel), and the reduced y 2 
for both treatments (lower panel) for our selected burst. We 
found that the f a increases to many times the pre-burst level, 
peaking at 17.8 ±4.7 at 0.25 seconds after the burst start. The 
errors are ler significance limits determined by XSPEC. The 
burst bolometric flux is calculated from the blackbody param- 
eters 



F = aT b \(R/d) 2 
= 1.076 x 10" 



kT bb V 
1 keV J 



K bb erg cm 2 s 



(2) 



where 7bb is the blackbody temperature, R is the effective ra- 
dius of the emitting surface, d is the distance to the source, and 
K bb is the normalization of the blackbody assuming isotropic 
emission. The calculated flux does not include the contribu- 
tion to the total flux due to the scaled persistent emission. The 
burst component flux for the variable f a fits is therefore con- 
sistently lower than the standard approach fits for the peak of 
the burst, where /„ > 1 . We also find that, after the beginning 
of the burst, the \ 2 for the variable f a fits is consistently lower 
than for the standard approach fit, with means of 1 .26 ± 0.59 
and 1.50 ± 0.60 respectively. A Kolmogorov-Smirnov test 
shows that the two sets of x 2 , values have a 4.0% probabil- 
ity of being consistent with one another. Figure 3 shows the 
flux-temperature curve for this burst. 

The pre-burst persistent emission is very much fainter than 
the peak of the burst emission, by a factor of 35. It is pos- 
sible that counting statistics in the burst spectrum could in- 
duce a spurious response in f a . To investigate this possibility 
we examined the spectrum with the highest f a , which was 
17.8 ±4.7 recorded 0.25s after the burst start. We took the 
blackbody parameters for that spectrum, IcT = 1.572 keV and 
normalization 1034 km 2 /(10 kpc) 2 , and generated 10 3 simu- 
lated spectra consisting of an absorbed blackbody with those 
parameters plus the persistent model times unity. The sim- 
ulated spectra also incorporated counting statistics typical of 
the detector. These represent a burst spectrum for which f a 
does not change and for which any measured deviation from 
f a = \ must be due to noise or the fitting procedure. 

We fit each of the simulated spectra with our variable f a 
model. The mean parameters for these were kT = 1.569 ± 
0.011 keV, K bh = 1061 ±49 km 2 /(10kpc) 2 and f a = 0.36 ± 
1.18. The mean parameters are therefore consistent with those 
that seeded the simulated fits. Adopting the standard deviation 
of the simulation values as the error on the fitted f a we found 
our measured f a for the real burst spectrum was 14.6 standard 
deviations higher than the simulated mean. We therefore con- 
clude that the high f a values do not arise from counting statis- 
tics alone. Figure 4 shows a histogram of the measured f a for 
the simulated spectra compared with the value measured from 
the real burst spectrum. 

Thus, for one burst we have found evidence for a spectral 
effect which can be described by varying the level of the per- 
sistent emission. This change is statistically significant for at 
least some of the spectra; 20 of the 37 of the spectra after 
the beginning of the burst had f a above unity with more than 
4er significance (using the above procedure, see also figure 4). 
The associated improvement in y 2 is of greater than 3tr sig- 
nificance using f-tests. We apply the new method to all the 
PRE bursts in our sample in the next section. 

4. STATISTICS OF MANY BURSTS 
4.1. Burst selection 

We restrict our sample of PRE bursts to exclude events that 
are unsuitable for analysis. Three bursts were excluded for 
which the Standard-2 data was missing, preventing estimation 
of the instrumental background via pcabackest. 

Some burst sources, such as 4U 1728-34, AQL X-l, 
and GRS 1747-312, lie in crowded fields containing other 
LMXBs. If the other source(s) were active at the time of 
observation, then their persistent emission could be confused 
with the source under observation and it would not be possible 
to scale only the persistent flux from the burst source. These 
sources need to be excluded from consideration. 
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FIG. 2. — Comparison of fitting a variable persistent emission factor, f a , to the standard approach fits for a burst from 4U 1636—536, where f a is a scaling factor 
of the pre-burst emission. The variable fit yields consistently lower burst fluxes (top panel), though these are within the error bars. The contribution to the flux 
from the variable persistent emission increases to many times the pre-burst level (middle panel) during the rise of the burst. The /„ factor remains elevated for 
a considerable time after the total flux begins to decrease. Allowing the persistent emission to vary improves the reduced \~ (lower panel). See Figure 3 for the 
changes in normalization and temperature for this burst. 
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FIG. 3. — Flux-temperature diagram for the 2001 Jun 15 burst from 
4U 1636—536, as measured by both the standard fit and the variable /„ fit. 
The fluxes shown are bolometric fluxes calculated from the blackbody pa- 
rameters, assuming a distance to the source of 6 kpc (G08). The variable fit is 
consistently slightly lower in flux, though within the uncertainties. Also plot- 
ted are lines of constant radius for 5 km, 10 km, 20 km, and 40 km. These 
radii do not take into account spectral corrections (e.g., Suleimanov et al. 
201 1). The radius of the photosphere at maximum expansion is the same us- 
ing both methods, but the variable f a fit gives a smaller, hotter photosphere at 
the touchdown point and beginning of the cooling tail. Typical error bars are 
given for points near peak flux; at lower flux they are much larger. 
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FIG. 4. — Histogram of the measured persistent emission factor f a for 1000 
simulated spectra for the burst from 4U 1636—536. These spectra consist of 
a blackbody with parameters taken from the real spectrum with the highest 
measured f a , plus the persistent model times unity. The scatter arises from 
the (simulated) counting statistics. The observed spectrum is 14.6 standard 
deviations from the predicted f a assuming the persistent emission does not 
change. 

We used All-Sky Monitor (ASM) data to quantify this ef- 
fect 9 . For each burst we found every source within the PCA's 
field of view and took the ASM counts at the time record near- 
est the burst start time. Since the PCA's response drops off 
approximately linearly with distance from the field of view's 
center, we multiplied the ASM counts by 1 - s, where s is the 
angular distance of the source from the center of the field of 
view, in degrees. We compared all the other sources to the 
burst source, and if the total contribution is more than 10% of 
the pre-burst flux of the burst source we excluded this burst 
from consideration. We found source confusion in 36 bursts. 

9 see xte.mit.edu 



The majority of these are from 4U 1728-34, which lies very 
near (A8 = 0.56°) in the sky to the Rapid Burster, and which 
was frequently observed to burst during targeted observations 
of the Rapid Burster (Fox et al. 2001). 

4.2. Modelling the persistent emission 

In order to measure a change in the level of persistent emis- 
sion for every burst spectrum we need to construct a model 
for the persistent emission for each burst that can then be in- 
corporated into the burst fits. It is not possible to simply use 
the detected photon counts for each energy bin because these 
include instrumental background. 

We found that the blackbody-plus-powerlaw model used in 
§3 does not give adequate fits for every burst, so we consid- 
ered a set of alternative spectral models. We fit each persistent 
spectrum with a range of different models. These are summa- 
rized in Table 1, and retained the fit that gave the best (ie. 
lowest) x\- We found that the six models listed in Table 1 
are sufficient to describe the persistent emission of all but one 
burst (see Figure 5). We are not overly concerned with theo- 
retical interpretations of these models- the main objective is 
to get a function that matches the data- but the literature does 
provide some commonly-used functions. White (1986) points 
out the importance of scattering, mentioning that this can be 
modelled with a Compton spectrum or a blackbody with an 
added power law, and we include both in our selection of al- 
ternatives. The inclusion of a Gaussian at 6.4 keV was moti- 
vated by the detection of fluorescent Fe emission for some 
sources at this energy (G08). The combination of thermal 
bremsstrahlung with a Gaussian was found by experimenting 
in XSpec with persistent spectra that could not be properly 
fit with any of the other five models. For models containing a 
Gaussian, the lower limit on the line width is set to 0. 1 keV to 
avoid the Gaussian simply approximating a delta function that 
removes the error on a single energy bin. In models where 
the hydrogen column densities are held constant, these are 
given in Appendix A. The Nh values are mostly taken from 
the literature and references are given in Appendix A. In per- 
sistent emission models where the hydrogen column densities 
are allowed to vary, the quantity A(E) (see equation 1) is to 
be thought of as as the product of the true interstellar absorp- 
tion and a multiplicative factor that corrects the shape of the 
persistent emission model. As A(E) is a multiplicative factor 
it is unaffected by changes in f a and so this treatment does 
not introduce any systematic effects. One persistent emission 
spectrum, preceding a burst from Cyg X-2 observed at 14:29 
March 27, 1996, (burst ID 2 in G08) could not be fit ade- 
quately by any of the six models, so this burst was excluded 
from the analysis. 

The reduced x 2 f° r the spectral fits to the remaining 332 
persistent spectra had a mean of 1.03 and a standard devi- 
ation of 0.35, with an average of v = 21 degrees of free- 
dom. The skewness was 1.00 and the kurtosis was 2.8. This 
distribution is therefore somewhat more skewed and signifi- 
cantly more peaked compared to expected theoretical values 
of \J$/v = 0.62 and Yljv = 0.57 for the skewness and kurto- 
sis respectively. However, a Kolmogorov-Smirnov test gave a 
value of D = 0.069 and a 8.8% probability that the measured 
distribution is consistent with the expected assuming a good 
fit. We therefore consider these persistent emission models 
adequate for use in the subsequent generation of burst fits. 
The distribution of xl for the fitted persistent emission spec- 
tra is shown in Figure 5. 
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FIG. 5. — Distribution of \\ for the fits to the persistent emission spectra for 
the 332 PRE bursts (blue stepped line) and the theoretical distribution (dashed 
curve) for v = 20. A Kolmogorov-Smirnov test showed a 8.8% probability 
that the recorded distribution is consistent with the expected assuming a good 
fit, indicating that our model fits of pre-burst spectra are adequate. 

4.3. Fitting burst spectra 

As in §3, initially we fit the burst spectra via the stan- 
dard approach using an absorbed blackbody, by subtracting 
the recorded pre-burst emission, ie. wabs*bbodyrad in 
XSpec. This is the same implementation of the standard ap- 
proach as the one outlined in G08. 

We then replace the detected pre-burst emission with just 
the instrumental background and fit the same spectra with 
a blackbody corrected by the adopted interstellar extinction 
value plus a multiple of the persistent emission model for that 
burst. We fix the hydrogen column density to the values given 
in Appendix A. The multiple of the persistent emission f a is 
allowed to vary as a free parameter, along with the tempera- 
ture and normalization of the blackbody, as we did in §3. 

The Levenberg-Marquardt algorithm used by XSPEC can 
become trapped in local % 2 minima, returning obviously un- 
physical results such as extremely high temperatures or nor- 
malizations. This problem can be mitigated if sensible initial 
values for the parameters are given rather than the XS PEC de- 
faults, which are not appropriate in every case. We use the 
blackbody temperature and normalization of the blackbody 
obtained from the standard model, and f a = 1, to seed these 
fits. 

Generating and fitting to 1,000 simulated spectra, as was 
done in §3.2, for each of the tens of thousands of burst 
spectra would be computationally prohibitive. To avoid 
needless work, we excluded spectra for which the vari- 
able f a fit was spurious. Only spectra with blackbody 
temperatures between 0.1 and 5.0 keV, normalizations of 
less than 10 6 km 2 /(10km) 2 , and bolometric fluxes less than 
10~ 6 erg cm" 2 s" 1 were retained. We also excluded f a determi- 
nations for which XSpec encountered fitting errors such as 
non-monotinicity or reaching the parameter limits. This left 
26,113 burst spectra out of the original 41,282. Almost all 
(>99%) of the discarded spectra were recorded before the be- 
ginning of the burst or very late in the cooling tail (i.e. stages 
or 4 in Table 2; see also §4.4), and the failure to obtain spec- 
tral fits can be attributed to very low photon counts for these 
spectra. We also reduced the number of simulated spectra for 
each measured spectrum from 1,000 to 320. 

As a further check that f a is measuring a real spectral effect, 
we took the highest f a spectrum for each burst and performed 
an f-test on it, comparing the x 2 statistic from the standard ap- 



FlG. 6. — Distribution of measured pre-burst persistent emission factor /„ 
for cooling tail spectra and Eddington limited spectra for all PRE bursts. The 
cooling tail spectra are strongly peaked around 1 , while the Eddington-limited 
spectra are not as strongly peaked. The distributions for the Eddington lim- 
ited spectra have more prominent tails, indicating a large number of elevated 
f a values, but there a significant number of high /„ values in the cooling tail 
also. The distribution for the expansion and cooling phases is consistent, sug- 
gesting that the fallback of the atmosphere does not contribute significantly 
to the emission. 

proach and variable f a fits. Because f a x P{E) is an additive 
component (see equation 1), f-tests are a suitable test (e.g., 
Orlandini et al. 2012). We found that, of these, 165 had de- 
tections of greater than 3a significance and 65 had detections 
of better than 4er significance. 

We define the burst stage: the pre-burst stage consists of all 
times before the beginning of the burst, the expansion stage 
occurs from the beginning of the burst to the time of maxi- 
mum normalization, the contraction stage is from the time of 
the maximum normalization up to and including the touch- 
down time, and the cooling tail stage is all times after the 
touchdown time but before the bolometric flux drops back be- 
low one quarter of the maximum flux. We neglect spectra af- 
ter this time. We refer to the expansion and contraction stages 
collectively as the Eddington-limited stage. These stages are 
summarized in Table 2. 

4.4. f a as a function of burst stage 

Figure 6 shows the distribution of persistent flux factor f a 
for every burst spectrum (ie. spectra from the Eddington- 
limited and cooling tail stages). The distribution of f„ has 
a wider spread for the Eddington-limited spectra, so that there 
is a larger fraction of high f a spectra during radius expansion 
than for the cooling tail. However, the population of high f, 
spectra in the cooling tail is not negligible. This suggests that 
elevated f a is not only a result of the burst being Eddington- 
limited. 

There are 45 1 spectra with f a less than zero, for which the 
f a measurement is significant to more than 3cr (via the proce- 
dure detailed in Figure 4), and which occur in the Eddington- 
limited phase, out of 26,113. Of these, 333 spectra come 
from just six bursts. These six events are summarized in Table 
3. By inspection of the flux-temperature curves we identified 
that four of the six are examples of superexpansion bursts (in't 
Zand & Weinberg 2010), particularly powerful PRE events 
for which the atmosphere is expanded to much larger radii 
than usual. In these events the temperature of the photosphere 
drops out of the detection band of RXTE, resulting in zero ob- 
served flux from the star. The other two are bursts with highly 
unusual flux-temperature curves and appear to consist of two 
consecutive expansions separated by a few seconds, followed 
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TABLE 1 

XSPEC MODELS FOR FITTING THE PERSISTENT EMISSION 



AjrEL mociei 


Number 


Notes 




of spectra 




wabs* (bbodyrad+power law) 


127 


nH fixed to values in Appendix A 


wabs* (bbodyrad+powerlaw+gauss) 


62 


nH fixed, Gaussian energy set to 6.4keV 


wabs* (compTT) a 


37 


nH fixed 


wabs* (bbodyrad+power law) 


18 


nH allowed to vary 


wabs* (gauss+bremss) 


31 


all parameters variable b 


wabs* (bbodyrad+diskbb c ) 


57 


all parameters variable b 


Total usable spectra 


332 




Rejected due to source confusion 


36 


Other active sources in field 


No background data 


3 




No good persistent model fit 


1 


Minimum \i > 5 



11 See Titarchuk (1994a) 

b All parameters are variable for the generation of persistent emission models. Their values are subse- 
quently frozen for the burst spectral fits in §4.3 

c See XSPEC manual (http://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/manual.html) and references 
therein 



TABLE 2 
Stages of a PRE burst 



Stage number Stage name Description 



Pre-burst Flux has not yet reached 1/4 of the maximum 

1 Expansion Normalization has not yet reached maximum 

2 Contraction Radius has reached maximum, kT has not yet reached maximum 

3 Cooling tail Flux still above 1/4 maximum 

4 Post-burst Flux has dropped below 1/4 maximum 



by a cooling tail with larger blackbody radius than the maxi- 
mum radius reached during expansion. 

4.5. Distribution of f a 

We compared the distributions of xt to the expected dis- 
tribution assuming a correct model using the Kolmogorov- 
Smirnov test. For the cooling tail we obtained D = 0.304 for 
the standard approach fit and D = 0.193 for the variable fit. 
For the Eddington-limited spectra we obtained D = 0.331 and 
D = 0.171 respectively. These results indicate that allowing 
f a to vary consistently improves the goodness-of-fit, and the 
improvement is more pronounced for the Eddington-limited 
spectra than for the cooling tail. These distributions are shown 
in Figure 7. However, none of the Kolmogorov-Smirnov tests 
were consistent with the null hypothesis; even our variable f a 
fits do not adequately fit the data, though they are an improve- 
ment on the old method. Some other systematic error must yet 
contribute to the discrepant spectral fits, perhaps deviations 
in the intrinsic burst spectrum from the assumed blackbody 
shape. Kolmogorov-Smirnov tests on the variable f a method 
against the standard approach likewise indicates that the two 
methods are not consistent with each other, for either the cool- 
ing tail or Eddington-limited spectra. 

We also investigated whether the degree of radius expan- 
sion correlates with the increase in persistent emission. We 
define the reduced radius for an individual burst as the photo - 
spheric radius divided by the radius at the touchdown point, 
which corrects for the distance of the source, assuming the 
isotropy factor does not change. We compared the maximum 
reduced radius with the maximum /„ for each burst; see Fig- 
ure 8. These two maxima do not necessarily occur at the same 
moment. We found only a slight relationship between the two 
quantities; a Kendall r rank correlation test found r = 0.138 
with 3.6cr significance. This suggests that obscuration of the 
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FIG. 7. — Distribution of reduced x f° r tne variable /„ vs. standard ap- 
proach spectral fits, for both cooling tail and Eddington-limited spectra. Also 
plotted is the theoretical distribution for 23 degrees of freedom, which is the 
mean for the standard approach fits (blue line). Allowing the persistent emis- 
sion to vary significantly improves the fits, but there are still significant de- 
viations from the expected distribution for a perfect model. The variable /„ 
and standard approaches are not consistent with being drawn from the same 
distribution using a Kolmogorov-Smirnov test. 
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TABLE 3 

Bursts with many negative/,, burst spectra 



Source 


Date and time 


G08 burst ID 


Description 


4U 1722-30 


1996, Nov 8, 07:00 


1 


Super-expansion burst" 


4U 2129+12 


2000, Sep 22, 13:47 


1 


Super-expansion burst 


2S 0918-549 


2008, Feb 8, 03:02 


5" 


Super-expansion burst 


4U 1722-30 


2008, Marl, 16:18 


4" 


Super-expansion burst" 


XB 1832-330 


1998, Nov 27, 05:45 


1 


Atypical burst profile 


HETEJ1900.1-2455 


2010, Sep 20, 05:29 


7 b 


Atypical burst profile 



a This burst was also studied by in't Zand & Weinberg (2010) 

b This burst postdates G08; a consistent burst ID numbering scheme is assumed 

c This burst was also studied by in't Zand et al. (201 1) 




Maximum burst r/r 

FIG. 8. — Maximum f a plotted against reduced radius (photospheric radius 
divided by radius at touchdown) for each burst. The two maxima do not 
necessarily occur at the same time. The dashed line shows the location of 
an accretion disc boundary layer at R = l.l&R* (Popham & Sunyaev 2001). 
The majority of photospheric expansion events exceed this radius. There is 
no obvious relationship between maximum f a and maximum expansion, nor 
is there a decrease in /„ outside the accretion disc radius, suggesting that 
obscuration of the disc by the expanding atmosphere does not significantly 
affect the observed accretion emission. 



inner parts of the disc does not greatly affect the observed 
persistent emission. 

One possible contribution to the elevated persistent flux 
during contraction is the fallback of the extended atmosphere. 
During the expansion phase, the atmosphere can be driven 
off at about 10 18 g s" 1 (Weinberg et al. 2006). This is, natu- 
rally, around the order of Medd, the Eddington accretion rate. 
As an order-of-magnitude estimate we make the assumption 
that half the expanded mass is driven off permanently. In a 
burst that is Eddington-limited for 5 seconds, the mass that is 
not expelled is about 2.5 x 10 18 g. There are approximately 
40% more contraction spectra than expansion spectra. All the 
spectra at these high fluxes are taken at 0.25s intervals, so the 
atmosphere must generally contract more slowly than it ex- 
panded, by a factor of about 40%. It follows that the accretion 
rate due to the atmosphere falling back could be as much as 
A^Edd/3. There are bursters for which the pre-burst accretion 
rate is only 1% of Medd, so for these sources there should be a 
noticeable excess of contraction stage spectra with f a around 
20-30. We have not detected any such excess (see Figure 6), 
so we conclude that atmosphere fallback does not contribute 
significantly to the accretion luminosity. A careful compari- 
son with non-PRE bursts would clarify this issue further, as by 
definition atmosphere fallback cannot occur in these events, 
and this will be investigated in a subsequent paper. 



5. SPECTRAL SHAPE CHANGES 

Deviations of the burst component (that is, everything ex- 
cept persistent emission) from a blackbody could also be re- 
flected by a changing f a , especially if these deviations man- 
ifest as an excess of photons at high energy. Such spectral 
changes have been theoretically predicted at both the high en- 
ergy (e.g., London et al. 1984, 1986), and low energy (e.g., 
Madej 1991; Madej et al. 2004) ends of the X-ray spectrum. 
Deviations from a pure black-body are frequently inferred 
during radius-expansion bursts (see, e.g., Kuulkers et al. 2002, 
2003). The influence of these deviations on f a measurements 
must be investigated, but this is difficult because the nature 
of the deviations is not known. Here we attempt to test the 
alternative hypothesis that the poor xt values given by the 
standard approach arise instead from deviations of the burst 
component from a blackbody. 

We take the approach of assuming the burst component re- 
tains a consistent spectral shape, which we model as a black- 
body plus power law. We returned to our "prototype" burst 
(see §3) and fit the Eddington-limited and cooling tail spectra 
(stages 1, 2, and 3 in table 2) with an absorbed blackbody plus 
powerlaw, with all parameters except Nh variable. The mean 
power law index of these fits was 1.93±0.68. We then fit the 
spectra again, this time holding the power law index fixed at 
this value but allowing the normalization to still vary. As we 
are here assuming that the power law is intrinsic to the burst 
emission, rather than being a separately varying component, 
we must specify the power law normalization so that it con- 
tributes a fixed proportion of the total flux. Using a linear fit, 
we found that the normalization of the power law was best 
described by 



Power law norm i 



Kbb x kT 4 
3212 



(3) 



where kT is the blackbody temperature in keV and is the 
normalization of the blackbody in km 2 /100kpc 2 . We fit the 
spectra a third time, this time tying the power law normaliza- 
tion to the blackbody parameters with the relation given by 
equation 3, which fixes the shape of the burst emission. 

We found the mean \i was 1.27 ±0.43 for the spectra of 
the Jun 15, 2001 burst from 4U 1636-536 using this fitting 
method, compared with 1.69±0.58 and 1.21 ±0.37 for the 
standard and variable persistent methods respectively. Thus, 
for this burst, the phenomenological model fits the data nearly 
as well as the variable persistent model. We then used the 
same model, with the same power law index and normaliza- 
tion relation, to fit every stage 1, 2, and 3 spectrum from all ra- 
dius expansion bursts in our catalog. Restricting these to spec- 



10 



tra where XS PEC did not encounter fitting errors for any of the 
three methods, we found that our phenomenological model 
gave a mean of 1.43±1.10, compared with 1.45±0.82 and 
1.21±0.60 for the standard and variable persistent methods 
respectively. The phenomenological model is thus globally no 
better than the standard approach. If there is a non-blackbody 
contribution to the burst component, it must differ between 
sources and change from burst to burst. Furthermore, it would 
also have to resemble the persistent emission for every burst, 
at least superficially, or the variable f, method would not be 
able to consistently improve the reduced x 2 ; and this itself 
hints at a relationship with the persistent emission. Clearly 
any further investigation into such deviations must be physi- 
cally motivated rather than phenomenological. 

We then tested whether introducing a phenomenological 
change to the burst spectrum removes the need for a variable 
persistent emission. We fit the spectra of the Jun 15, 2001 
burst from 4U 1636-536 using the variable persistent model, 
but replacing the blackbody burst component with the phe- 
nomenological spectrum. If our initial detection of a variable 
f a merely reflected the non-blackbody character of the burst 
component, then we would expect the detection to disappear. 
We still detect /„ to vary significantly, though the values are 
consistently lower by roughly 20% than for the original fits. 
Furthermore, the mean xt f° r these fits was 1.17±0.35, sug- 
gesting that allowing the persistent emission to vary improves 
the fits even when a non-blackbody model for the burst com- 
ponent is assumed. These results suggest very strongly that 
our detection of f a cannot be attributed to a confounding spec- 
tral effect intrinsic to the burst component. 

If our variable f a fits are measuring a non-blackbody con- 
tribution to the burst spectrum then we should expect f a to 
increase with burst flux. Suppose that a non-Planckian part 
of the burst component contributes a fixed fraction of the 
burst component flux and that our variable f a fits are trying 
to remove it. Then if the burst flux doubles, we would re- 
quire /„ to double also in order to fit out the non-Planckian 
part. We investigated this by comparing f a against F/FEdd, 
where F E dd is the source's Eddington flux (see Appendix A) 
for all Eddington-limited and cooling tail spectra. It is clear 
from Figure 9 that there is no relationship between the quan- 
tities. To test whether plotting all the points together ob- 
scures a relationship present in individual bursts we selected 
the Eddington-limited and cooling tail spectra from our pro- 
totype burst and three randomly selected other bursts, and 
performed Kendall rank correlation tests on their f a against 
their blackbody fluxes. No statistically significant correlation 
was present in any of these bursts. If the hypothesized non- 
blackbody contribution is present at the high energy tail, then 
we would expect an inverse correlation between blackbody 
temperature and f a because at low temperatures the blackbody 
drops out of the band detectable by RXTE, leaving only the 
hard tail to be fit. However, we find a slight positive corre- 
lation of r = 0.031 at 4.2cr significance using a Kendall rank 
correlation test. 

Our new fitting model holds the shape of the persistent 
emission fixed. Since the persistent emission of LMXBs 
changes shape according to whether it is in a high or low state 
(Hasinger & van der Klis 1989), it is reasonable to think it 
may change shape temporarily if the accretion rate changes 
during a burst. It is possible that spectral shape changes in the 
persistent spectrum contribute to the high \ 2 m our approach. 
However, it is not clear that an accretion rate enhancement 
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FIG. 9. — f a against normalized burst flux for all Eddington-limited and 
cooling tail spectra. There is no obvious correlation between the two quanti- 
ties. There is a clearly visible clump at low flux and f a ~ 1, showing that f a 
returns to its original level after the burst. There are also many points with 
high fa and normalized fluxes between 0.2 and unity, indicating that f a can 
remain elevated a considerable way down the cooling tail. 

due to radiation drag will have the same effects on the per- 
sistent spectrum as the neutron star's usual movement around 
the Z-track, since they have very different physical causes. 
Indeed, in't Zand et al. (2013) found in their study of a burst 
from SAX J1808. 4-3658 that the persistent spectrum became 
harder, whereas ordinarily they would expect increased accre- 
tion to soften the spectrum. 

6. COMPARISON WITH THEORY 

W92 performed ID simulations of geometrically thin, ax- 
isymmetric irradiated accretion discs around a neutron star. 
Their models include the effects of viscosity and general rela- 
tivity. W92's results predict that radiation torque from bursts 
can enhance the accretion rate by up to two orders of mag- 
nitude. If our f a is assumed to be entirely an accretion en- 
hancement, then it corresponds identically to their quantity 
AM* /M* and is predicted to be inversely correlated with pre- 
burst accretion rate, accretion disc viscosity parameter j3, neu- 
tron star spin frequency, and neutron star radius. All of W92 's 
models assume a neutron star of mass 1 .4M Q . 

Table 1 of W92 lists peak accretion enhancements for the 
computed models; in Figure 10 we compare the maximum f, 
measured for each burst against 7, the pre-burst accretion flux 
as a fraction of the Eddington flux of the burst source. The 
Eddington flux for each neutron star is the mean peak flux 
of every PRE burst observed from that star; see Appendix A 
for the values and details of their calculation. Following G08 
we measure 7 by integrating the chosen pre-burst persistent 
model between 2.5 and 25 keV and dividing the resulting flux 
by the burst source's Eddington flux. W92's Table 1 lists re- 
lated quantities. W92's models begin with a non-rotating neu- 
tron star with radius 9km, 7 of roughly 0.3, and disc viscosity 
parameter, /?, of 10~ 4 . They then allow the spin frequency, 
accretion rate, j3, and radius to vary in turn, while holding the 
other quantities fixed at their original values. 

The disc viscosity parameter f3 (Coroniti 1981) is similar 
to the Shakura-Sunyaev disc viscosity (Shakura & Sunyaev 
1973) but relates the viscosity to the gas pressure rather than 
the total pressure, which differs in discs which are radiation- 
pressure dominated. W92 gives accretion rates in dimension- 
less units: m = Mc 2 /LEdd, where M is the mass accretion rate 
and Ledd is the Eddington luminosity, whereas we give accre- 
tion rates in terms of the energy release. Since the mass M* 
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and radius of the neutron stars in W92 are specified, we 
have 



Our results in Figure 10 show a decrease in peak f a with 
increasing 7, and the slope is consistent with the predictions 
of W92's three points representing a nonrotating star with in- 
creasing accretion rate, while our peak /„ values are signif- 
icantly lower than those predicted for a nonrotating neutron 
star. If we assume a moderate (5 and rotation frequency of 
about 300 to 600 Hz (e.g., Muno et al. 2001; see also Ap- 
pendix A), then not only the observed correlation of peak f a 
with 7, but also the normalization, appear to be consistent 
with W92's predictions (see Figure 10). The models predict 
that stellar rotation period and disc viscosity parameter /3 have 
a large influence on peak f a but that neutron star radius appar- 
ently has little influence. The upper edge of our measured 
points appears to be very roughly consistent with jf a < 1, 
implying that the accretion luminosity cannot greatly exceed 
Lgdd- This is consistent with the predictions of Burger & Katz 
(1983) and Miller & Lamb (1996), who find that M Edd is a 
natural cap on the accretion rate. 

In Figure 1 1 we show maximum f a against spin rate for all 
neutron stars in our sample whose spin rate is known. These 
are listed in Appendix A. W92's models predict a gradual de- 
crease of maximum f a with spin frequency. Our data shows 
maximum f a values consistent in magnitude to W92's mod- 
els, but it is difficult to discern any trend, because of the large 
scatter in individual sources (due to the additional dependence 
on 7 and the disc viscosity) and the small range of known neu- 
tron star spin frequencies (since these are all accreting neutron 
stars). We would require PRE bursts from slowly rotating 
neutron stars to better constrain this relationship. A further 
difficulty is that W92's models do not reach the Eddington 
limit, so we are comparing PRE bursts in observations to non- 
PRE bursts in models. 

As the W92 models are one-dimensional, care must be 
taken in applying them to accreting neutron star systems in 
which three-dimensional effects are likely to be important. 
Further theoretical studies, preferably in three dimensions, 
would be very valuable. 

7. DISCUSSION 

We have performed an observational investigation into 
whether or not the persistent emission contribution varies dur- 
ing photospheric radius expansion bursts. We allowed the pre- 
burst emission to vary, parametrized by a factor which we de- 
note /„. We detected a statistically significant increase in f a 
for nearly all the PRE bursts in the catalog, suggesting an en- 
hanced (rather than suppressed) level of persistent emission 
during a burst. Our new method results in a significant im- 
provement in the reduced \ 2 of spectral fits compared to the 
standard model but not to the level of formal statistical con- 
sistency for all the spectra. 

Since the persistent emission is known to be an approxi- 
mate measure of the accretion rate onto the neutron star, we 
interpret our results to indicate that the accretion rate onto the 
neutron star generally increases during bursts. Obviously the 
persistent spectrum could also change shape in response to 
a varying accretion rate, and there are suggestions (Homan 
et al. 2007) that the X-ray luminosity may not be proportional 
to the mass accretion rate. There are, however, currently no 
predictions about how the persistent emission spectrum may 
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FIG. 10. — Maximum f a against 7 (the pre-burst accretion rate as a fraction 
of Medd) for all PRE bursts (black squares). Also plotted are the results of ten 
computer simulations performed by W92 (see their Table 1). The real bursts 
show generally lower f a , but the same slope, as W92's nonrotating neutron 
star model. If a moderate fi and spin frequency of 300 to 600 Hz is assumed, 
then our results appear to agree with W92. The shaded region indicates the 
approximate area spanned by W92's models. The bottom edge of the shaded 
box is uncertain due to the fact that none of W92's models vary both rotation 
rate and disc viscosity simultaneously, but it may extend further down than 
shown here. 
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FIG. 1 1 . — Maximum f a against neutron star spin frequency for all PRE 
bursts arising from sources whose spin frequency is known (see Appendix 
A). Also plotted are the results of three computer simulations performed by 
W92 with spin frequencies of 10~ 5 , 333, and 667 Hz. Although there is large 
scatter due to the dependence of /„ on pre-burst m and p, our results are 
consistent in magnitude with W92's predictions. The anti-correlation of peak 
f a with spin period predicted by theory is not obviously visible in our data. 

change in response to rapidly increasing radiation drag, so we 
make the simplest assumptions that it does not change shape 
and that f a oc M. We have shown that peak f a measured dur- 
ing each burst anti-correlates with the pre-burst accretion rate 
with a slope consistent with that predicted by the theoretical 
models of Walker (1992), who investigated the effect of radi- 
ation drag on the accretion disc. If the effects of neutron star 
spin are accounted for, the magnitudes of the observed peak 
f a are also consistent with theory. This suggests that our de- 
tection of an increased pre-burst persistent emission reflects 
an enhanced accretion rate due to radiation torques on the ac- 
cretion disc around the star. 

All of W92's models are sub-Eddington, so some care must 
be taken in extrapolating their results to PRE bursts. However, 
W92 points out that the expanding atmosphere in a PRE burst 
should have very little angular momentum compared with the 
accretion disc. Radiation coupling between disc and envelope 
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can thereby spin down the disc, increasing the accretion rate. 
This would also have the effect of spinning up the envelope 
and boosting its expansion. Such a mechanism may help ex- 
plain why the atmosphere takes longer to return to the neutron 
star surface than to reach maximum radius. 

It may be argued that a deviation of the burst spectrum from 
a blackbody could mimic a variable accretion rate by intro- 
ducing a high energy excess that our method then attempts 
to remove. If the deviation is of fixed size then any spu- 
rious f a measurement it causes will also anti-correlate with 
the pre-burst accretion rate 7, similarly to Figure 10. How- 
ever, we would then also expect f a to be constant with con- 
stant burst flux, but we found it to vary greatly during the 
Eddington-limited phase when the burst flux is approximately 
constant. We found no anti-correlation between f a and black- 
body temperature, which would occur if the blackbody com- 
ponent drops out of RXTE's detection band and leaving only 
the deviation to be fit out. We also attempted to model a hard 
tail deviation in the burst spectrum using a power law with 
fixed index and normalization tied to the burst flux. This did 
not improve the spectral fits for the entire collection of bursts, 
and did not cause our detection of f a to be suppressed. While 
in general we do not expect to see discrete spectral features in 
PRE bursts, in superexpansion bursts they can be visible (in't 
Zand & Weinberg 2010). We have found that in these events 
f a consistently drops rather than rises, consistent with zero 
flux from the source. As suggested by in't Zand & Weinberg 
(2010) this may be due to the emission from the burst com- 
ponent dropping out of the band detectable by RXTE together 
with the persistent emission being obscured by the superex- 
panding shell, though it is also possible that the vigorously 
expanding envelope disrupts the accretion disc and thereby 
temporarily halts accretion. An enhanced f a cannot be at- 
tributed to the presence of spectral features superimposed on 
the burst component continuum. While we have attempted 
to exclude confounding spectral effects intrinsic to the burst 
component as an explanation for f, enhancements, we cannot 
rule out the possibility of other interpretations. 

If our variable persistent fits causes the inferred photo- 
spheric radius at touchdown to differ systematically from the 
standard fits, then this would have implications for determi- 
nations of the equation of state of neutron star matter, such 
as those of Ozel et al. (2009); Giiver et al. (2010); Steiner 
et al. (2010). We investigated this by taking the ratio of the 
touchdown radii as determined by the variable f a and standard 
methods for all the PRE bursts in our sample and obtained a 
mean value of 0.97 ±0.11. This indicates that the inferred 
touchdown radius can differ by ~ 10%, and that there is a 
slight trend towards lower neutron star radii using the vari- 
able persistent method. We have found that f a frequently re- 
mains elevated after the end of the Eddington-limited phase; 
this may have implications for studies that use the cooling 
tail to constrain the neutron star parameters (e.g., Galloway & 
Lampe 2012). 

7.1. Structure of the accretion disc 

The details of the transfer of material from the inner edge 
of the accretion disc to the neutron star surface are still uncer- 
tain (e.g., Bildsten 1998). It is not known if it occurs at a close 
inner boundary layer, or if there is some mechanism that regu- 
lates the infall. PRE bursts offer a means of settling this ques- 
tion: to be subject to radiation pressure the expanding photo- 
sphere must be optically thick and therefore must obscure our 
view of everything interior to it. The atmosphere expanding 
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FIG. 12.— A non-PRE burst from 4U 1636-536 observed on August 10 
2005. The f a factor is significantly elevated, suggesting that the level of 
persistent flux rises for non-PRE bursts as well as PRE bursts. The time 
evolution of /„ is similar to the one shown in 2, indicating that the physical 
process is similar in both cases. 

so far that it covers all of the emitting parts of the accretion 
disc would cause the persistent emission to be reprocessed 
inside the optically thick atmosphere and effectively become 
part of the burst emission. Even very modest radius expansion 
would hide the boundary layer due to this "shrouding", and 
cause fa to decrease, since the boundary layer is thought to 
be geometrically small and located at a radius ~ 1.2/?* (e.g., 
Popham & Sunyaev 2001). Thus our finding of consistently 
high f a during a burst argues against the existence of a thin 
boundary layer that remains near the surface of the star and 
dominates the persistent emission (see Figure 8), but is con- 
sistent with a boundary layer that becomes much wider in re- 
sponse to increased luminosity as suggested by Popham & 
Sunyaev (2001). 

We also note that performing the analysis method of this pa- 
per on non-PRE bursts may provide further information about 
the structure of the accretion disc. In these bursts the atmo- 
sphere does not expand and will not obscure any of the accre- 
tion disc. A careful comparison of non-PRE bursts against 
PRE bursts, which progressively obscure parts of the disc, 
might therefore give insights into structure and properties of 
the disc. A preliminary investigation of a number of non- 
PRE bursts using the same analysis method indicates that the 
f a values increase during these events as well (see figure 12), 
which would be expected if increased M is caused by radiation 
drag as in W92's models, but which would not be expected if 
f a is tracking spectral changes caused by the emitting photo- 
sphere distending. A proper study of non-PRE bursts will be 
the subject of a subsequent paper. 
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8. APPENDIX A 

Table 4 shows the Nh values we adopt for each burst source 
studied in this work, the reference from which we obtained 
that value, and the number of photospheric radius expansion 
bursts recorded by RXTE from that source, including bursts 
we discarded due to source confusion or unsuitable spectral 
data (see §4.2). 

We also give Eddington fluxes, measured using the standard 
approach, for all sources that have PRE bursts recorded by 
either RXTE (i.e. in the extended G08 catalog) or the Wide 
Field Camera on BeppoSAX 11 . The Eddington fluxes are the 
means of the peak fluxes for every PRE burst recorded for 
that source, weighted by the inverse square of the error of the 
measurement. 

If only one PRE burst has been recorded for any source, we 
simply report the peak flux for that burst and the error on that 



measurement. For sources with two or more PRE bursts we 
calculate a reduced Xmeas assuming the source has constant 
Eddington flux. If Xmeas i s consistent with a constant flux, 
we allow the flux to vary up and down such that it now has 
Xnew = 1 + Xmeas' tne amount by which we vary it is the error 
on the original measurement. If x meas is not consistent with 
constant Eddington flux we artificially scale the errors on the 
individual measurements until x meas = 1 . Then we allow the 
flux to vary such that Xn ew = 2- 

The prolific burster 4U 1636-536 has PRE bursts that are 
known to be bimodal in flux (Ebisuzaki & Nakamura 1988; 
Galloway et al. 2006). Most of its PRE bursts are thought to 
take place in a pure helium atmosphere and have peak fluxes 
around (68.6 ± 5.5)x 10~ 9 erg cm -2 s -1 . A few bursts have 
peak fluxes lower by a factor of m 1 .7, and these are believed 
to take place in a hydrogen-rich atmosphere (G08). We con- 
sider both regimes separately. The lower flux bursts include 
a tentatively identified PRE event recorded by RXTE on 2000 
Jan 22, 04:43:48 UT. We include this burst to calculate the 
mean Eddington flux for these low flux events, but exclude it 
for the rest of the analyses in this paper. 

We also list neutron star spin frequencies, where known. 
These have been calculated either from burst oscillations (BO) 
(e.g., Watts 2012), or from pulsar timing. 
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TABLE 4 

Nh values chosen for spectral fits to burst data, Eddington fluxes, and spin periods 



Source 




N H (10 22 cnT 2 ) 


No. of PRE bursts 


Ref. 


Eddington flux 


Spin (Hz) 


Method c 


Ref. 






RXTE (WFCf 


(Nh) 


10~ 9 erg cm -2 s _1 






(Spin) 


/itt ncn , 1 r\ 

4U 051 J— 40 


0.03 


4 (3) 


1 


14.5 ± 3.5 








hXO 0748—676 


0.80 


5 (0) 


2 


46.5 ± 4.6 


552 


BO 


27, 28 


2S 09 1 8-549 


0.35 


2(3) 


3 


1 19.2 ± 12.4 








4U 1608—522 


0.89 


17 (5) 


4 


172.2 ± 21.8 


620 


BO 


27, 29, 30 


4U 1636-536 (nigh) 


0.25 


76 (2) 


5 


7.1.9 ± 6.8 


58 1 


BO 


27, 31, 32 


4U 1636-536 (low) 


0.25 


2 b (0) 


5 


41.5 zfc 1.2 


581 


BO 


27, 31, 32 


MXB 1658-298 


0.20 


12 (0) 


6 


17.0±4.0 


567 


BO 


27, 33 


XTE 1701-462 


2.00 


2 (0) 


7 


43.4 ± 1.4 








4U 1702-429 


1.87 


5 (3) 


8 


87.7 ±4.5 


329 


BO 


27, 34 


4U 1705-44 


1.90 


4 (0) 


9 


41.0±3.8 








XTE J1710-281 


0.40 


3 (0) 


10 


7.1 ± 1.5 








4U 1722-30 


0.78 


3 (23) 


1 


61.7 ± 12.4 








4U 1728-34 


2.60 


94 (0) 


1 1 


95.0 ±8.4 


363 


BO 


27, 35 


KS 1731-260 


1.30 


4(3) 


12 


48.6 ±5.6 


524 


BO 


27, 36, 37 


4U 1735-444 


0.14 


10 (0) 


13 


34.2 ±5.6 








GRS 1741.9-2853 


11.30 


6(1) 


14 


35.3 ± 10.9 


589 


BO 


27, 38 


1A 1742-294 


1.16 


2 (0) 


8 


37.8 ± 1.4 








SLX 1744-300 


4.50 


4(3) 


15 


13.9±3.1 








GX 3+1 


1.59 


2 (0) 


16 


60.0 ± 1.4 








SAX J1748.9-2021 


0.47 


11 (0) 


1 


38.0 ±6.0 


442 


Pulsar 


27, 39 


EXO 1745-248 


3.80 


2(0) 


1 


69.0 ±2.8 








4U 1746-37 


0.26 


3 (0) 


1 


5.4 ±0.7 








SAX J1747.0-2853 


8.80 


11 (2) 


17 


52.5±7.1 








IGR J 17473-2721 


3.80 


3 (0) 


18 


1 13.5 ± 12.1 








IGR J 17498-2921 


1.20 


1 (0) 


8 


51.6±1.6 


401 


Pulsar 


27, 40 


XTE J 1759-220 


2.84 


3 (0) 


8 


15.7±0.8 








SAX J1750.8-2900 


0.90 


2(1) 


8 


54.1 ±2.1 


601 


BO 


27, 30, 41 


GRS 1747-312 


1.39 


3 (0) 


1 


13.4±4.4 








oAX J loUo.4— 365o 


0. 12 


8 (2) 


19 


230.1 ± 13.2 


40 1 


Pulsar 


27, 42 


XTEJ1810-189 


4.20 


1(0) 


20 


54.2 ±1.8 








SAXJ1810.8-2609 


0.35 


1(1) 


21 


111.2±3.0 








GX 17+2 


1.90 


2(0) 


22 


15.5 ±0.5 








4U 1820-303 


0.16 


16 (34) 


1 


60.5 ±4.0 








XB 1832-330 


0.05 


1 (0) 


1 


33.7 ±4.4 








HETEJ1900.1-2455 


0.16 


7(0) 


23 


123.9 ±10.6 


378 


Pulsar 


27, 43 


Aql X-l 


0.40 


14(0) 


24 


99.6±21.3 


550 


Pulsar 


27,44 


XB 1916-053 


0.32 


12 (0) 


25 


30.6 ±3.6 


270 


BO 


27,45 


4U 2129+12 


0.03 


1(1) 


1 


40.8 ±1.6 








Cyg X-2 


0.05 


8 (0) 


26 


13.1 ±2.1 








Ser X-l 


0.38 


6(0) 


26 


29.4 ±7.1 
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